Spin Dependent Joule Heating due to Rashba Coupling and Zitterbewegung 
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Investigating microwave absorption in asymmetric Si quantum wells in an external magnetic 
field, we discover a spin dependent component of Joule heating at spin resonance. We explain this 
effect in terms of Rashba spin-orbit coupling which results in a current induced spin precession and 
Zitterbewegung. Evidence is based on the observation of a specific dependence of the electron spin 
resonance line shape and its amplitude on the experimental geometry which in some range suggests 
a "negative" differential power absorption. 
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I. INTRODUCTION 

In solids with lower than mirror symmetry, the 
Rashba-type of spin-orbit coupling for charge carri- 
ers in semiconductors and metalaDZI gives rise, e.g., to 
Dyakonov-Perel spin relaxation, the spin Hall effect, elec- 
tric dipole spin resonance or a group of spin photo- voltaic 
effects.!^ In this paper we present evidence for another 
effect originating from Rashba coupling. We show that 
it opens a new channel of spin dependent energy trans- 
fer between an r.f. electric field and the kinetic energy 
reservoir. 

Rashba coupling has recently received attention also 
since, making use of spin-orbit (S-0) interaction, efficient 
spin manipulation may be achieved by exciting spin pre- 
cession by an r.f. electric field. This mechanism can be 
much more efficieniP (one of the main preconditions for 
spintronic^^Hl) than the classical excitation by an oscil- 
lating magnetic field which is used in standard electron 
spin resonance, ESR. Also, addressing nanostructures via 
electric currents is much easier than the generation of an 
r.f. magnetic field in a small volume. 

Spin excitation via the S-0 field appears well 
understood.'^ For low frequency, the excitation efficiency 
is determined by the electric current, i.e., it scales with 
the electron mobility (Current Induced ESR).^" For high 
frequency, when the displacement current exceeds the 
drift component, momentum dissipation can be neglected 
and spin-fiip transitions can be treated within the formal- 
ism of electric dipole transitions. Within this limit, this 
effect has been described many years ago by Rashba and 
it is called Electric Dipole ESR (ED ESR)." 

The magnitude of the S-0 field has been well docu- 
mented in Si quantum wells-fi^ There is, however, a rich 
and hardly understood dependence of the CI ESR line 
shape and its anisotropy on exp erime ntal condition^^ 
in Si and also in other materials.l^^lSl in particular, all 
models based on the effective S-0 field predict that the 
ESR line shape can be described by the imaginary-, i.e., 
the so-called absorption component of the dynamic mag- 
netic susceptibility, x {^) — both for the S-0 compo- 



nent and for a ny superposition of S-0 — and magnetic 
r.f. fields .12^ This conclusion is in contrast with ob- 
servation: the power absorption spectra are composed 
of absorption- and dispersion-like componentJ^ of xi^)- 
Moreover, the amplitude of the dispersive component 
(proportional to $He(x) = X can be much bigger 

than the absorptive one. This effect is puzzling also be- 
cause X i'-^) changes its sign in some frequency range 
which seems to suggest the possibility of a "negative 
power absorption" . 

In this paper, we investigate the ESR line shape for 
both electric and magnetic microwave excitation of 2D 
electrons in SiGe/Si/SiGe quantum wells and we propose 
a theoretical explanation based on the fact that the S-0 
field results in a dependence of the electron velocity on 
spin orientation. Consequently, excitation of spin pre- 
cession leads to oscillating components of the electron 
velocity. The interference of this motion with the r.f. 
electric field modifies the electric absorption in the sam- 
ple contributing thus to the power absorption spectrum. 
This effect is proportional to the precession amplitude 
and thus it appears at ESR. This specific type of electri- 
cal absorption can be treated thus as arising from a spin 
dependent electric conductivity or as a resonant modifi- 
cation of the Joule heat. In agreement with the exper- 
imental conditions, we model the effect within the lin- 
ear response limit assuming small precession excitation, 
where saturation effects in the spin polarization can be 
neglected. In particular, we omit higher order effects, 
such as a dependence of the electric conductivity on spin 
polarization, i.e., the so-called polarization signal.^ 



II. EXPERIMENTAL 

We investigate conduction electron ESR in a one- 
sided, modulation doped Si quantum well defined by 
Sio.75Geo.25 barriers as described in Refs. 5, 12. The 
same type of samples, grown by molecular beam epitaxy 
on relaxed Sio.75Geo.25 layers, is used here. We use a 
Bruker ElexSys E580 ESR spectrometer in the X-band at 



2 



temperatures of 2-40 K. Two types of microwave cavities 
were employed: the standard one (TE201) where sam- 
ple is placed in the maximum of the magnetic microwave 
field, Hi, and a special cavity (TE301) where the sample 
is located in a maximum of the electric field, Ei. In both 
cases El is parallel- and Hi perpendicular to the exter- 
nal static magnetic field, H. The sample can be rotated 
about the direction of Hi , which allows to change the di- 
rection of H with respect to the sample growth direction, 
n, described an angle 6. The microwave bridge works in 
such a way that the observed signal monitors the power 
absorption spectrum only. Because of the modulation of 
the external field, the spectrum is obtained as the first 
derivative of the power absorption. 



III. EXPERIMENTAL RESULTS 

The ESR due to the 2D conduction electrons in Si 
produces a very narrow line (typically << 1 G) super- 
imposed on a broad spectrum of electrical absorption. 
The latter is caused by the excitation of cou pled plasmon- 
cyclotron oscillations and harmonics .li^liSI Here the spec- 



tra are rather complex, showing up to 4 resonances^ in 
the accessible range of 0-1.5 T. The width of these reso- 
nances allows to estimate the momentum relaxation time. 
Results spread in the range of Tp — {3- 10~^^-4 • 10"-^^) s. 

A typical example for the ESR line shape is shown 
in the inset of Fig. [l] It is very far from the classical 
absorption line (dashed line) but it can be well fitted by 
a sum f{H) = A f (H) + A f {H) of imaginary and 
real parts derived from the Lorentzian shape function: 

/±H = ±-^^^^^^ = ./±H + i/±M- (1) 

Here T2 is the transverse spin relaxation time and 
wl = 7-ffESR the Larmor frequency at ESR, and 7 is the 
gyromagnetic factor. For the ESR line shown in the inset 
of Fig. Ill the amplitude of the dispersion component, A , 
is by a factor 2.3 bigger than A . Both amplitudes, and 
their ratio change from sample to sample. They depend 
also on the direction of the applied field, 9, as shown in 
Figs. [1] and [2| 

The data in Fig. [l] were obtained with the sample po- 
sitioned in the node of Ei. The classical ESR absorption 
(induced by Hi ) is well evidenced only for small by a 
dominant absorption amplitude, A (see Fig. [l]). Here, 
close to its node, Ei is small and roughly perpendicu- 
lar to the sample layer. Therefore almost no current is 
induced. Consequently all contributions caused by the 
electric field are small, and the absorption signal origi- 
nates from magnetic dipole (MD) transitions only. 

The MD ESR signal is expected to be independent of 
the sample orientation. But when 9 is increased, an in- 
plane electric current is induced. Accordingly, A should 
increase, but it decreases. Obviously another signal ap- 
pears, compensating part of the amplitude of MD ab- 
sorption. Simultaneously, the amplitude of the dispersive 
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FIG. 1: Angular dependence of the amplitudes of the real, A 
(red dots), and the imaginary part, A (blue squares), of the 
ESR signal at = 27r- 9.4 GHz of a high mobility 2D electron 
gas {uTp = 3) when the sample is placed at the node of 
El. Inset: observed derivative of microwave absorption ESR 
line (solid green line). The line shape is fitted (dashed green 
line) by a combination of real, / (B) (dash-dotted red), and 
imaginary, / (B) (dotted blue), Lorentzians (Eq. [l]). In this 
example, A is by a factor 2.3 bigger than A . 



component, A , substantially increases. Both compo- 
nents scale roughly like {Ei sin 9)^, indicating that both 
are caused by the S-0 field.fi^ 

Placing the sample in the node of Hi, we obtain the 
angular dependencies of the signal amplitudes shown in 
Fig. |2] In that case, for — >■ both amplitudes practi- 
cally vanish and increase with 9. Obviously, in that case 
the amplitude of the MD resonance is negligible and all 
observed ESR signals are excited by Ei sinO via S-0 cou- 
pling. The inset shows the derivative of CR absorption 
spectra, i.e., the dynamic electric conductivity, measured 
for 9 = 0°,30°,60° and 90°. The width of CR for that 
sample indicates a comparatively low mobility. With in- 
creasing 9, the CR moves towards higher field, roughly 
with l/cos0, demonstrating the dependence of CR on 
the perpendicular component of the static magnetic field, 
H cos 9. 



IV. MODELING AND DISCUSSION 

The occurrence of a dispersive component in the ab- 
sorption spectra (which sometimes is even negative) can- 
not be explained by existing models,'^'^ neither for MD 
ESR nor ED ESR, nor for CI ESR. But the experimen- 
tal data clearly indicate that the effect is induced by the 
r.f. electric field, demonstrating that it is related to S-0 
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coupling. 

We consider here a system lacking mirror symmetry. 
In that case, a zero-field spin splitting occurs, described 



by the Rashba^^^term, % — aYi{sxky — Syk^), or equiva- 
lently, by the S-0 field, , acting on individual electrons 
with a momentum Hk: 



• s = aR(k X n) • s. 



(2) 



where s stands for its spin, or is the Bychkov-Rashba co- 
efficient, g the g-factor, fiB the Bohr magneton and fxo is 
the vacuum permeability. In the presence of the Rashba 
term in the Hamiltonian the velocity of each electron is 
modified: 

V fc =;^ = ^^,r =— + ^ nxs , 3 
at n m n 

The second term on the r.h.s. of Eq. ^ depends on 
spin. A precessing electron thus undergoes an oscillatory 
motion, which has recently been discussed as an example 
for Zitterbewegung.^^^^ 

The electron current j is obtained by averaging the 
spin-dependent velocity over the k-space, yielding: 
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(n X M). 



(4) 



The microwave drift (or displacement-) current, jd(i^) = 
(7Ei(cj), described by the first term, in the low field limit 
is proportional to the microwave electric field Ei(i) = 
Die [Ei(w)e^"^*] and the complex dynamic conductivity, 
a{uj,H). The second term contains the magnetization, 
M = ns5/iB (s). Here we omit its static component Mq 
because we investigate the microwave power absorption 
resulting only from spin precession. The latter causes 
oscillating components, and thus the second term also 
contributes to an oscillating current. Within the linear 
response approximation the amplitude of the precession 
is given by Af±(a;) = X±{'-^)H±{itj), with the dynamic 
magnetic susceptibility x±{'-^) — ""T-^o/iCw). Here the 
indices ± stand for the amplitudes in a rotating coordi- 
nate system and thus: M±(w) = (Ma(w) =F ^/^{uj)) 
where a, /3 designate axes perpendicular to the direction 
of H. The amplitude of the oscillating field H± (oj) , which 
drives the precession is composed of the microwave mag- 
netic field, Hi±{Ld), and the S-0 field, -ffR±. The mean 
value of the Rashba field is obtained by averaging Eq. (d) 
for all electron spins 
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(5) 

The first term results from the displacement of the Fermi 
sphere caused by a drift current for ujTp << 1 (valid also 
for a displacement current for uTp » 1). The second 
one is a consequence of a change in the spin subband pop- 
ulation which causes a non-zero mean k-vector in each 
spin subband even without a drift current. Consequently 
this term contributes to the Rashba field and it is re- 
sponsible e.g. for the spin Hall effect.'^ However, a fast 
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FIG. 2: Angular dependence of the amplitude of dispersive 
(dots: experimental, circles: evaluated from CR spectra) and 
absorption (squares, experimental) components for a sample 
of a moderate mobility {uuTp — 0.3), placed in the node of 
Hi. Inset: first derivative of the CR absorption spectra for 
different angles 0. 



precession {ujT2 >> I) docs not change the occupation 
of spin states. Therefore, the precessing magnetization 
does not result in any additional Rashba field and only 
the static (or some low frequency component) of magne- 
tization, Ms, may occur in Eq. ([5]). Summarizing, the 
oscillating amplitude H^± (w) is determined by the oscil- 
lating current amplitude, jd±(w), and described by the 
first term of Eq. 

Finally, the spin related, oscillating component of the 
total current (second term in Eq. Q) scales with a^, 
the dynamic susceptibility, X±('^)i Jd±('^)- There- 
fore, both contributions to the current are proportional to 
Jd±('^) and thus to E±{io). But they differ in phase. The 
drift current jd(i) is shifted in phase with respect to the 
electric field Ei{t) by (pjE = arctan (wr^), where Tp is the 
momentum relaxation time. The magnetization induced 
current is additionally shifted by a difference in phase be- 
tween Unit) and M(t). According to the solution of the 
Bloch equations this shift is ifHM — arctan ((w — WL)^2)■ 
The electric power, i.e., the Joule heat, absorbed per 
unit area F of the electron gas is given by dP/dF — 
(j(i) • Ei(t))j, where the current density is obtained from 
Eq. Q. Because of the interplay of these phase shifts 
the spin dependent contribution is a linear combination 
of imaginary and real part of the magnetic susceptibility, 
and the amplitude of both components are described by 
various components of the electric conductivity tensor. 

Putting the a;-axis into the sample plane, the in-plane 
component of the microwave electric field is Ex{oj) = 
Eiiuj) sin0. Then the spin dependent contribution to the 
Joule heat is: 
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dF 



C 



{[91, 



(6) 



r 



with C = a^m* / fiog"^ n'^h'^ns. Here we neglect contri- 
butions related to the ESR inactive component, x_(a;). 
Eq. ([6| is the main result of this paper. It shows a spin 
dependent contribution to the Joule heat with a reso- 
nance at the Larmor frequency. This type of signal has 
both an absorptive and a dispersive component. The ex- 
pressions in the square brackets are proportional to the 
two amplitudes, Aj and Aj , respectively, of the ESR sig- 
nals originating from the spin dependent motion of the 
electrons. Both amplitudes can be positive or negative, 
depending on the electron conductivity. A negative sign 
shows only that the total Joule heat at resonance is par- 
tially reduced by the spin dependent component. 

Because of the complex behavior of the electric con- 
ductivity of a high mobility 2D electron gaJ^^ modeling is 
hardly possible, but for low mobility samples (wTp << 1), 
the first term, Aj, can be approximated by the first term 
in the square bracket in Eq. ^ . In that case, the ampli- 
tude Aj is expected to be proportional to the CR absorp- 



tion at ESR, i.e., to dP{ujo)/dF = ^Rz{(t^^{ujo)) \E^\ /2. 
Therefore, the angular dependence of A (8) can be di- 
rectly compared with the experimentally measured dy- 
namic conductivity at ESR, ^^^{d , H^gji) , which is ob- 
tained when the sample is placed in a node of Hi. 
The open circles in Fig. [2] represent the product of 
sin^ 8 ■ ^Rccr^xiO, H^sr)- This dependence fits the experi- 
mentally observed angular dependence of A (8) very well, 
proving the proposed origin of the dispersive component 
in the absorption spectra. 



A similar analysis of A {8) is much more difficult. 
There are two contributions due to the energy exchange 
between the microwave and (i) the kinetic energy of the 
electrons and (ii) their magnetic energy: the former, Aj , 
is described by the second term in Eq. ^ and the latter 
by the amplitude of CI ESR, A-^, which is proportional 
toP 



dPM(w) _ C 



dF 



4 e^n, 



[Icr^^P + 2cos83m{a*yaxx) + cos^ 8\axy\^] wx+H • EI{uj). 



(7) 



The transfer to the magnetic energy (ii) has only an 
absorption component, proportional to x ■ The exper- 
imentally accessible amplitude of the absorption signal 
is obtained by the sum: A — A^^ + Aj . Neither of 
them can be derived from the measured electric absorp- 
tion without detailed modeling of the absorpt ion j'i^ which 
is beyond the scope of this paper. 

Therefore we discuss only a qualitative picture by ap- 
proximating the electric conductivity within the Drude 
model. For that case we show the expected angular de- 
pendencies of the signal amplitudes in Fig. [3] for wTp — 1. 
The peculiarities seen close to 8 — 80° correspond to the 
coincidence of the spin- and cyclotron resonances for an 
effective electron mass of m* = 0.2mo. They become 
more pronounced for bigger values of WTp . 

Within the Drude model the amplitude of the disper- 
sive signal is expected to have a positive sign only (see 
Fig. k5|). The absorption components A^ and A^ are 
of opposite sign and their sum is expected to change 
sign at some direction of the applied field. Altogether a 
quite complicated angular dependence results except for 
61 90°, where Aj'(90°) = -^^^(90°), i.e., the absor 



tion signals are expected to compensate each other.E^ 



All amplitudes resulting from the Drude model are 
expected to scale with uJTp. The inset in Fig. [3] shows 
the dependence of the ESR signal amplitudes on ujTp 
for 6 = 90°. Within the limit of drift conductivity (for 
LUTp << 1) the absorption amplitudes are proportional to 
the square of electron mobility and to (ojTp)^ whereas Aj 
, is proportional to WTp. Therefore Aj dominates the ab- 
sorption components. In the opposite limit of a displace- 
ment current, the absorption amplitudes become inde- 
pendent of frequency and momentum dissipation. There- 
fore, in the infrared range, the absorption signal is ex- 
pected to exceed Aj which decreases with (ajrp)~^. In 
fact, the spectra observed in the high frequen cy ra nge 
usually exhibit an absorption type of line shape-lS^Ell 



V. CONCLUSION 

In conclusion, we presented a new efi^ect resulting from 
the Rasba type of S-0 coupling which manifests itself 
by a new type of spin resonance. The r.f. Rashba field, 
caused by an r.f. current, induces spin precession, leading 
to an increase of the magnetic energy. There is, however. 
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FIG. 3: Angular dependence of the different ESR contribu- 
tions according to Eqs. ([6||7| assuming the Drude model for 
the conductivity. The inset shows their dependence for in- 
plane orientation of H on uiTp. Here the static magnetization 
is assumed to be proportional to frequency. 

also a component of the electric current which depends on 
spin orientation. Therefore the induced spin precession 
leads to an additional, oscillating component of the cur- 



rent, modifying the Joule heating. This spin dependent 
current is shifted in phase as compared to the classical 
drift current, and thus the Joule heating can increase or 
decrease at spin resonance, depending on the interplay of 
phases. The spectral dependence of the magnetic channel 
is given by the absorption shape function and predicts a 
positive absorption only, but the line shape of the elec- 
trical channel is given by a combination of absorptive 
and dispersive components where the amplitude of the 
absorption usually has a negative sign. Consequently, 
the absorption signals of the two channels partially com- 
pensate each other and the resulting amplitude can be 
of either sign. Anyway, the absorption spectrum with 
the dispersive component is an attribute of the electric 
channel, evidencing Zitterbewegung at the Larmor fre- 
quency. Our model predicts thus a unique dependence of 
the resonance spectrum and of the amplitude on the ex- 
perimental geometry. The experiments confirm this well, 
in spite of the fact, that details of the dynamic electric 
conductivity of small samples of high mobility 2D elec- 
tron gas cannot be easily modeled. 
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